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ABSTRACT 


Fifteen lidar observations of the stratospheric aerosol were made 
between February and November 1975* All observations revealed the 
greatly increased particulate backseat ter ing that followed the eruption 
of the volcano Fuego in October 1974. Vertical structure consisted ini- 
tially of multiple layers, which later merged to form a single, broader 
peak. Essentially all of the increased scattering was confined to alti- 
tudes below 20 km. Hence, aerosol layer centroi-is in 1975 were typically 
several km below their altitude prior to the eruption. 

Our observations began in mid-February, at about the time of maximum 
northern midlatitude influence of the volcanic injection. From late Feb- 
ruary on. both vertically integrated particulate backseat ter ing and the 
peak ratio of particulate to gaseous backscattering displayed approxi- 
mately exponential declines, with mean lifetimes (l/e-lives) of eight 
and eleven months, respectively. These relatively short residence times 
are a combined consequence of the low altitude of the volcanic particles 
and their larger mean size as compared to the preinjection, or unperturbed, 
aerosol. The peak scattering ratio of our average 1975 profile was 1.7, 

and the vertically Integrated particulate backscattering was 3.6 x 1C 
-1 

sr (both at X ■ 694 nra). The mean midvlsible particulate optical thick- 
ness, derived from measured backscattering and realistic optical models, 

li 


was about 0.03, approximately six times the mean value in the year before 
the Fuego eruption, but not as large as values observed for some years 
after the 1963 Agung eruption. 

Radiative and thermal consequences of the measured post-Fuego layer 
were computed using several recently published models. The models pre- 
dict a temperature increase of several °K at the altitude of the layer, 
caused by the infrared absorption bands of the sulfuric acid particles. 
The surface temperature decrease predicted by the models is considerably 
smaller than 1°K, partly because of the small optical thickness of the 
volcanic layer, and partly because of its short residence time relative 
to the earth-ocean thermal response time. 
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1. Introduce ion 

In October 1974 the volcano Fuego in Guatemala (14.^*N, 91° W) 
erupted violently over a period of several days (Smithsonian, 1974). 

The eruption plume penetrated into the stratosphere, and very shortly 
thereafter dramatic twilight effects were observed at many northern 
hemisphere locations as the cloud was transported by the stratospheric 
circulation (e.g., Meinel and Meinel, 1974; Fegley and Ellis, 1975; 

Volz, 1975). Sudden increases in the numbers, mass, and optical back- 
scattering and extinction of stratospheric particles were measured by 
balloon-borne counters, aircraft -borne samplers, ground-based Ildars, 
and a spacecraft -borne radiometer at a variety of locations (e.g., Hof- 
mann and Rosen, 1976; Ferry and Lem, 1975; Lazrus et al. , 1975; Remsberg 
and Northam, 1975; McCormick and Fuller, 1975; Fegley and Ellis, 1975; 
Fernald and Frush, 1975a, b; Hlrono et al., 1975; Pepin, 1975). At 
northern midlatitudes, both the numbers of particles and their optical 
backscattering were observed to increase very irregularly until about 
February 1975, after which a more regular decline began. 

This paper describes stratospheric lidar measurements made at 
Menlo Park, California (37° N, 122? W) between mid-February and mid- 
November 1975. The measurements thus began when the influence of the 
Fuego injection was near its peak in northern midlatitudes, and they 
trace the subsequent decline and tendency toward pre-Fuego conditions. 
The measurements were an extension of a pre-Fuego observation series 
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that has been described elsewhere (Russell et al., 1974, 1976a, b). 

The purpose of this paper is to document the 1975 decline in stratospheric 
aerosol content, to compare our post-Fuego lidar measurements to pre- and 
post-Fuego measurements made by a variety of methods, and to derive some 
conclusions about the radiative and thermal consequences of this strato- 
spheric aerosol event. Section 2 reviews the equipment and data analysis 
procedures used in the study. Section 3 presents the results and compares 
them to other measurements. Section 4 discusses the significance of the 
results, especially in terms or particle transport and removal processes, 
and possible radiative /climatic consequences of the Fuego or similar 
events. 
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2. Instrumentation and Data Analysis 

A lidar operates by firing a pulse of laser energy into the atmo- 
sphere and measuring the amount of energy scattered back by atmospheric 
constituents. The ground-based ruby lidar used for all observations in 
this study is the one used in our pre-Fuego series and was described in 
detail by Russell et al. (1974; 1976a, b). It produces 1- to 2-Joule 
pulses of wavelength 694 nm and uses a 40.6-cm (16-inch) receiving tele- 
scope. The receiver is interfaced to a computer-based digital data ac- 
quisition system that permits real-time display of signals integrated 
from many pulses. For each nightly observation, the lidar was operated 
in four or more alternating periods of low and high sensitivity, with 
the receiver partially obscured to produce the low-sensitivity mode. 

This method was used to extend the vertical range over which received 
signals could be accurately digitized. Both current and pulse-counting 
data were acquired to extend the vertical range of digitization accuracy 
within each sensitivity mode. 

The final data analysis was performed off-line on a larger computer 
us*ng the method described by Russell et al. (1974; 1976a, b). Data from 
all four operating modes ( 'ov-sensivity current, low-sensitivity pulse- 
counting, high-sensitivity current, and high-sensitivity pulse-counting) 
were combined to obtain a signal profile valid over an altitude range 
that typically extended between ~7.5 and 30 km (and above). This sig- 
nal profile, together with a molecular density profile measured by 
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radiosonde at Oakland (32 km to the north) on the night of the Ildar 
observation, and a stratospheric extinction profile, was used to compute 
a profile of scattering ratio, defined as 

f (a) + t (z) 

R(z) « ■ * " T V T T i ♦ * (*)/f <*> . (1) 

£ \Z) p in 

in 

where f (z) and f (z) are, respectively, the particulate and molecular 
p m 

(gaseous) backscatterlng coefficient 0 at altitude z. 

Because the Ildar system calibration and the low-altitude atmospheric 

transmission are unknown, the scattering ratio profiles can initially be 

computed only to within an arbitrary constant. To fix this constant, 

each profile of R(z) is normalized to produce a statistically significant 

minimum value of unity. This procedure yields valid profiles, provided 

a level of negligibly small particulate-to-molecular backscatterlng ratio 

(f /f ) exists within the altitude range of signal analysis (see Section 4 
p m 

for discussion). In computing scattering ratio profiles, a particulate 
extinction profile must be assumed, and in this work the model of Slterman 
(1968) for a wavelength of 700 nm was used. This model differs from the 
model of negligible particulate extinction used in our 1972-74 observa- 
tions, and is used to t'*ke account of the increased stratospheric partic- 
ulate content observed since the eruption of Fuego (see also Section 4). 
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The probable error of each computed ecattering ratio data point ie 
also computed ae a part of the analyst*. As described in detail by 
Russell et al. (1974; 1976a, b) the computed error includes the probable 
error in lidar-* sasured signal, in radio sonde -measured density, and in 
assumed stratospheric extinction. "Error bars," equal to ± one standard 
error derived In this manner, are then plotted along with the measurement 
data, as shown below. 
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3. Results 


The vertical profiles of scattering ratio measured in the fifteen 
1975 observations are shown in Figure 1* Also shown for comparison in 
each frame Is the average profile of scattering ratio measured between 
June 1973 and March 1974 In our pre-Fuogo series. In each profile the 
excess of scattering ratio over unity may be Interpreted as a type of 
"optical mixing ratio," as it is equal to the ratio of light backscattered 
by particle* to that backscattered by gases [cf. Eq (1)]. Thus, for ex- 
ample, on 13 February at 19 km, particulate backscatterlng was three times 
an great as gaseous backscatterlng for the 694 nm wavelength. This very 
thin and intense scattering layer was evidently quite localized or short- 
lived, since a scattering ratio as large as four was not measured in sub- 
sequent observations either by us or by other North American observers. 
Nevertheless, the major increase in stratospheric particulate content 
caused by the Fuego eruption is obvious in all of the 1975 observations. 

It is evident that the injection had its largest effect at altitudes 
below 20 km, thus shifting the peak altitude of scattering ratio at least 
several km below its immediate ore -Fuego value of 22 to 23 km. Other 
features of aerosol layer vertical structure revealed by the lidar mea- 
surements Include a multiple-layer structure (in the early 1975 obser- 
vations) that gave way ,:o a single, broader peak; a temporary increase in 
particulate backscatterlng at the lowest stratospheric altitudes in 
spring; and a change in layer-cop structure (near 20 km) from very abrupt 
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FIGURE t VERTICAL PROFILES OF SCATTERING RATIO 
MEASURED BY THE RUBY LIDAR IN 15 
OBSERVATIONS. FEBRUARY-NOVEMBER 1975 

The dashed profile, shown for comparison, is the 
average of 16 observations made in the pre-Fuego 
period, June 1973-March 1974. 
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(i.e. "flat") to more gradual, or sloping, and then back to somewhat 
abrupt in November 1975. All of these features are in agreement with 
the balloon observations made in Laramie, Wyoming, during the same period 
by Hofmann and Rosen (1976). 

In the later 1975 profiles there is a definite increase of scattering 
ratio above 24 km, suggesting a gradual penetration of volcanic influence 
(either particles or particle* forming gases) to these altitudes. It is 
especially interesting to note that the November profiles, and to a cer- 
tain extent the 3 September profile, have a distinct minimum in scattering 
ratio above the main peak. Above this minimum a very broad secondary 

layer can be seen, as shown in more detail, and at higher altitudes, 
in Figure 2. These profile shapes, which cannot be caused by the normal- 
ization procedure, suggest a particle formation source that is above the 
main layer. Recently, Crutzen (1976) has described photodissociation of 
(possibly volcanic) carbonyl sulfide (CSO) as a possible source of sulfur 
above 25 km, which would lead to particle formation at those altitudes. 
Although the present data sample is very limited, its crude agreement 
with Crutzen' s theory warrants further examination (see also Section 4). 
[Scattering ratios above 35 km are shown as dotted lines in Figure 2 
because local radiosonde data on the observation dates did not extend 
above that altitude. At altitudes between 35 and 40 km, possible differ- 
ences between actual and standard atmospheric densities are large enough 
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FIGURE ? VERTICAL PROFILES OF SCATTERING RATIO MEASURED 
BETWEEN 11 AND 13 NOVEMBER 1975 

Dotted curves, in the region above the radiosonde density 
measurement, were computed from the midlatitude spring-faM 
standard atmosphere 
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(~157.) that scattering ratios computed from a standard density profile 
are not useful for normalization or for the purposes of Figure 2.} 

One reason for making the four November observations in a period 
of only three days was to determine the size of day-to-day changes in 
layer structure. As can readily be seen from p igures 1 and 2, such 
changes were small, and caused primarily by vertical movements of the 
layer as a whole, indicating that layer structure had by November become 
rather stable. Within the November data set, dawn and twilight observa- 
tions were made to permit a search for significant structure differences 
occurring between these two times of day. Again, close comparisons did 
not reveal any dawn-to-twilight or midnight-to-twilight differences. 
Finally, on an even shorter time scale, four individual portions of the 
13 November dawn run, each consisting of several tens of minutes of data 
selected from a two-hour period, were compared in a search for systematic 
differences that might occur as sunlight began to impinge on the aerosol 
layer. This search also proved negative, as any differences that were 
observed did not appear to have any systematic dependence on the increas- 
ing solar illumination. 

Figure 3 shows the time variation during 1975 of three different 
quantities that measure the predominant altitude of particulate back- 
scattering between 10 and 30 km. These quantities are the altitude of 
maximum scattering ratio; the centroid of the "optical mixing ratio," 

R-l; and the centroid of the particulate backscattering coefficient f^. 
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FIGURE 3 TIME VARIATION OF ALTITUOE OF MAXIMUM 
SCATTERING RATIO. CENTROID OF (R-1) AND 
CENTROID OF PARTICULATE BACKSCATTERING. 
FEBRUARY-NOVEMBER 1975 

The centroid calculation! were restricted to the altitude 
range 10-30 km. The altitude of R mM in the June 
1973-March 1974 average profile is 21.9 km, or 0.5 km 
above the dashed line. 


As expected from Figure 1, all three altitudes are consistently lower 
than the corresponding values observed during 1973-74. The springtime 
increase of scattering at the lowest stratospheric altitudes, mentioned 
above, is reflected by a lowering of centroids between late February and 
late April. After mid-April there is a gradual overall rise in centroids, 
as they appear to be approaching the 1973-74 values. Such a rise in cen- 
troids can be attributed to a reduction of scattering at lower altitudes 
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compared to higher altitudes , and is consistent with the occurrence of 
more efficient particle removal mechanisms at lower altitudes than at 
higher. As discussed by Hofmann, et al. (1975) and Hofmann and Rosen 
(1976), there are two possible mechanisms that tend to produce this result. 
One mechanism is the variations in tropopause height between late spring 
and early fall, which provide an effective cleansing process that is not 
present above 17 km. The second mechanism is particle sedimentation. 

This mechanism was more effective at lower altitudes than at higher, be- 
cause in the post-Fuego period, as shown by Hofmann and Rosen (1976), 
the mean particle size between 14 and 18 km was larger than it was at 
higher altitudes. Faster sedimentation of larger particles therefore 
leads to the observed result. Again, the presence of a particle- 
production mechanism above the initial injection, as suggested by Crutzen 
(1976), could also contribute to this effect. 

Figure 4 shows the behavior of maximum scattering ratio and verti- 
cally integrated particulate backscattering (10 to 30 km) over the 

<j|( 

February -November 1975 period. An overall decline is evident in both 


The values for 26 February and 30 April have been increased somewhat 
(relative to the values in Figures 1 and 2) by a renormalization that 
forces scattering ratios above 22 km to approximate those measured in 
surrounding observations. This procedure is to compensate for the 
original normalization above 10 km, a region that the other lidar ob- 
servations indicate frequently does not contain a "clean" layer suitable 
for normalization. 
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FIGURE 4 TIME VARIATION OF VERTICALLY INTEGRATED 
PARTICULATE BACKSCATTERING AND MAXIMUM 
SCATTERING RATIO. FEBRUARY-NOVEMBER 1975 

The vertical integration was restricted to the altitude 
range 10-30 km. 


quantities. After the first February observation, the maximum scattering 
ratio values (open circles) are fitted fairly well by an equation of the 
form 


R - 1 = 
max 


( 


H? 

ma v 



-t ft 
e 


( 2 ) 
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where Che best-fit mean life, t, is ll months. This value converts to 
a half-life t^ of 7.7 months. Prior to April, the vertically integrated 
particulate backscattering values (dark circles) do not decrease as regu- 
larly as do the (R • l) values (open circles). Nevertheless, the mean 
max 

life of an exponential best- fit provides a useful measure of the overall 
rate of decline. The values thus obtained are a mean life t of 8.4 months 
and a half-life t^ 2 of 5.8 months. Thus vertically integrated particu- 
late backscattering declined more rapidly than (R -l). Again, this 

max 

is consistent with more effective removal mechanisms below the altitude 

of R , because the majority of vertically integrated particulate back- 
max 

scattering occurred at those altitudes. (Note again the centroid alti- 
tudes in Figure 4 and the discussion thereof.) 

The relation of the present postvolcanic maximum scattering ratios 
to those measured by a number of observers since 1964 is shown in Fig- 
ure 5a. The February and 9 April 1975 values of (R -l) were a factor 

max 

of eight or more larger than typical values observed at SRI between June 
1973 and March 1974, shortly before the Fuego volcanic eruption. Never- 
theless, typical values since March 1975 have not been as large as those 
frequently observed by Grams and Fiocco (1967) in the years following 
the Agung eruption of March 1963 (see also Grams 1966). Since the 
18-month observation series of Grams and Fiocco did not even begin until 
nine months after the Agung eruption, these data indicate that less 
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1962 1964 1966 1968 1970 1972 1974 


SA-4019-12R1 


FIGURE 5 COMPARISON OF HOAR ANO BALLOON-BORNE 
PARTICLE COUNTER MEASUREMENTS OF THE 
STRATOSPHERIC AEROSOL. 1963-1975 

a. Maximum ratio of particulate to gaseous backscattering 
measured by a number of stratospheric lidar groups. Boxes 
indicate mean of many observations; bars indicate range in 
time and magnitude. GF: Grams and Fiocco (19671, 
Massachusetts; CL: Collis and Ligda (1966), California; 

C. Clemesha et al. (1966), Jamaica; S: Schuster (1970), 
Colorado; CR: Clemesha and Rodrigues (1971), Brazil; 

0: Ottway (1972), Jamaica: FS: Frush and Schuster 
(unpublished), Coloiado; F: Fox et al. (1973), Hawaii and 
Bermuda; YE: Young and Elford (1975), Australia; RVH: 
Present work and Russell et al. (1974, 1975b), California. 

b. Number of particles (radius ^0.15 pm) above tropopause 
as measured by photoelectric particle counter (Hofmann 

et at.; 1972-74; 1973b. 1974, 1975, 1976). Arrows give 
times of volcanic erruptions with appreciable stratospheric 
penetration. 
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particulate and particle- forming material was Introduced into the northern 
mid-latitude stratosphere by Fuego than by Agung. 

Figure 5b shows values of the number of particles above the tropo- 
pauae, as measured by balloon-borne counters flown from several Northern 
midlatitude sites (see references in the Figure caption). The general 
trend of these data parallels that of the lidar data, in showing an over- 
all decline by an order of magnitude between 1964 and late 1974, an ab- 
rupt increase between October 1974 and early 1975, and a quite regular 
decline thereafter. Considering the differences between the two quanti- 
ties being plotted ( a point mixing ratio in Figure 5a and a vertically 
integrated absolute concentration in Figure 5b) and also the differences 
in measurement location, the overall agreement between the two data sets 
is very good. 

It would of course be desirable to compare measured quantities 
that are more similar in character. However, until recently, lidar data 
have not been published as vertically integrated backscattering amounts, 
and balloon data have not been summarized as maximum mixing ratios. For- 
tunately, this situation is changing, and Hofmann and Rosen (1976) have 
recently published a time series of the maximum mixing ratio (particles 
per mg of air) and vertically Integrated content (15 to 20 km), both as 
measured by balloon, in a format very similar to that of Figure 4. The 
agreement between the lidar and balloon data revealed in this manner is 
quite striking. Both the balloon-measured maximum mixing ratio and the 
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balloon-measured column content displayed approximately exponential de- 
clines after mid-February 1975, as did the lidar data of Figure 4. More- 
over, the maximum mixing ratio had a longer mean life (12 months) than 
did the column content (6.6 months), again as in Figure 4. It may be 
significant that the balloon -measured mixing ratio lifetime was somewhat 
longer (by 8%) than the corresponding lidar-measured value (11 months). 

The difference could be caused by the stronger weighting of Ildar back- 
scattering by larger particles, which tend to fall out faster than small 
ones. The Ildar and balloon column-content measurements cannot be com- 
pared exactly, because of their different height ranges. 

Figure 6 shows the error-weighted average profile of scattering 
ratio for the February-November 1975 period. (The 26 February and 30 
April profiles were omitted from the average because their normalization, 
above 10 km, makes them prone to underestimation. Moreover, only a single 
average profile from November was included, to avoid biasing Che 1975 
mean toward year-end values.) The dotted curves show the typical obset- 
vational uncertainty in a single nightly observation. The dashed curves, 
which show the standard deviation of the profiles Included in the average, 
indicate the large degree of variability during the nine-month period of 
observation, and show that most of this variability was confined to 
altitudes below 20 km. 
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FIGURE 6 AVERAGE VERTICAL PROFILE OF SCATTERING 
RATIO MEASURED OVER CENTRAL CALIFORNIA 
COAST FFBRUARY-NOVEMBER 1975 

Solid line is error-weighed mean of 13 observations 
Dotted lines indicate typical experimental uncertainty 
in a single profile Dashed lines indicate standard 
deviation of observations included in data set 


The February-Novembor 1975 average profile of scattering ratio Is 
compared to our average 1973-early 1974 profile and the 1964-65 Crams 
md Fiocco (1967) profile in Figure 7. The large increase of particulate 
backscatter Ing In the recent post-Fucgo period relative to the intnedlate 
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SCATTERING RATIO 

SA-4019-10R1 

FIGURE 7 AVERAGE VERTICAL PROFILES OF SCATTERING 

RATIO OBSERVED IN THREE RUBY LIOAR STUDIES 
OF THE STRATOSPHERIC AEROSOL BETWEEN 
1964 AND 1975 

For references see Figure 5. 


pre-Fuego period Is readily evident, as Is the predominantly lower alti- 
tude of the post-Fuego scatter!**.?. Nevertheless, the even larger average 
scattering and lower predominant altitude of the post-Agung observations 
are equally evident. 

A characteristic of the stratospheric aerosol layer that Is of pri- 
mary Importance to climatic studies Is Its optical thickness, or vertically 
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integrated extinction# Derivation of an optical thickness from the present 
Ildar backscattering results requires the use of a particulate backscatter- 
to-extlnction ratio (backscattering phase function), which has never been 
measured for the stratospheric aerosol. Nevertheless, useful estimates 
of optical thickness can be derived by using model backscattering phase 
functions computed for realistic optical models of the stratospheric aero- 
sol. Recently, several such models, based on measurements, have been 
published (Russell et al., 1976a, b; Cadle and Grams, 1975, Plnnlck et 
al«, 1976; Toon and Pollack, 1976). Most models agree on a particle 
composition of concentrated sulfuric acid, except possibly just after 
volcanic injections. Since recent measurements Indicate that the post- 
Fuego aerosol also was sulphuric acid (Ferry and Lem, 1975; Hofmann and 
Rosen, 1976), we adopt this composition in the discussion to follow. 

There is less agreement on particle size distributions, probably 
because of different measurement methods and because of actual variability 
in this characteristic. (A post-Fuego Increase in mean particle size has 
been well documented by Ferry and Lem (1976) and Hofmann and Rosen (1976). 
Possible size distributions Include the Haze H and Haze L models of 
Deirmendjian (1969, 1973), a lognormal distribution derived by Plnnlck 
et al (1976) from balloon-borne optical particle counter measurements, 
and a zeroth-order logarithmic distribution (Zold) derived by Toon and 
Pollack (1976) from a rather extensive set of measurements. Backscatter- 
ing phase functions for these size distributions (for a 700-nra wavelength 
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and a 757. sulfuric acid refractive Index of 1.42- 01) are, respectively, 
0.013, 0.018, 0.015, and 0.017 sr~* (a range of ± 16% of the mean). 

The range of turbidity profiles that there values yield when com- 
bined with our single 1975 average backscattering profile (Figures 6 and 
7) is indicated by the shaded area in Figure 8. [Turbidity is here de- 
fined as the ratio of particulate to molecular extinction. The profiles 
have been converted to a wavelength X of 550 nm assuming that particulate 
extinction varies as X \ the approximate dependence displayed by the 
optical models of Pinnick et al. (1976).] Also shown for comparison is 
the turbidity profile derived by Elterman (1975) from searchlight measure- 
ments made in New Mexico in November 1974 (during the imnediate post-Fuego 
increase) by using an assumed bistatic scattering phase function (Reeger 
and Seidentopf, 1946). The magnitudes of the lidar and searchlight pro- 
files are similar, as might be expected from the time dependence shown 
in Figure 5. However, the altitude of the searchlight-measured turbidity 
maximum is about 3 km below the lidar maximum. This also is attributable 
to the temporal development of the Fuego injection, and is consistent with 
both lidar and balloon observations made in November and December 1974 
(Fernald and Frush, 1975; Hofmann and Rosen, 1976). 

Table 1 lists values of particulate optical thickness derived from 
the present and previous measurement series. As can be seen, the mean 
particulate optical thickness during the 1975 post-Fuego decline was 
about 8 ix times as large as the value in the immediate pre -Fuego years. 
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TURBIDITY 

SA-4019-20 


FIGURE 8 PROFILES OF STRATOSPHERIC TURBIDITY (RATIO OF 
PARTICULATE TO GASEOUS EXTINCTION) DERIVED 
FROM THE PRESENT SERIES OF LIDAR OBSERVATIONS 
ANO FROM A SEARCHLIGHT MEASUREMENT IN LATE 
1974 (ELTERMAN, 1975) 

The shaded region is derived from the single average scattering 
ratio profile shown in Figures 6 and 7, plus a range of model 
backscattering phase functions. 
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Table I 


PART1CUIATK OPTICAL THICKNESS, t, BElfcBBN 10 
FOR A WAVELENGTH . - S30 i 


KM AND LISTED ALJITOOK, 



0.38-0.52 


Elterraan (1968). Based on average of blstatic searchlight measurements and model phase 
function. 

Y 

Plnnlck et al. (1976). Based on numerous balloon-borne particle counter measure- 
ments and several realistic optical models. 

* Russell et al. (1976a, b). Based on average of 16 Ildar measurements and a range of 
model backscattering phase functions. 

Toon and Pollack (1976). Based on tidar and sky-brightness measurements and model 
phase functions, plus trend of balloon, searchlight, and solar transmission 
measu rcments. 

* 

Elterman (1975). Based on searchlight measurement and model phase function. 


Present work. Based on average of Ildar measurements and range of model backscattering 
phase functions. 

















4. 


Discussion 


Temporal development of aerosol structure 

The temporal development of the northern midlatitude post-Fue go 
stratospheric aerosol observed in this study is the result of several 
interacting processes: transport of particles and gases from the orig- 

inal injection point (14. ^N, 91°W), vertical diffusion from the injec- 
tion strata, gas-to-particle conversion, particle growth and sedimenta- 
tion, and loss of particles below the tropopause. When this study began, 
in February 1975, the northern midlatitude stratospheric aerosol concen- 
tration was near its postvolcanic maximum, evidently as the combined 
result of northward transport of the initial injection cloud and of 
gas-to-particle conversion processes. A multiple-layered structure 
was evident in most of the February to April 1975 observations, evidently 
a lingering result of the original injection structure noted by many 
observers (e.g. , Hofmann and Rosen, 1976; Hirono et al. , 1975; Fernald 
and Frush, 1975; McCormick et al.,1975). By May the multiple layers had 
merged to form a single, broader layer, the centroid of which rose with 
time toward the altitude of the preinjection centroid. As discussed in 
Section 3, this upward centroid movement is consistent with the occurrence 
of faster particle removal at lower altitudes than at higher, as well as 
upward diffusion from the injection strata. The faster removal of 
particles at low altitudes (10-17 km) could be caused by vertical tropo- 
pause movements (Hofmann et al., 1975, Fernald and Frush, 1975b) or by 
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faster sedimentation of the larger particles which occurred preferen- 
tially at those altitudes (Hofmann and Rosen, 1976). 

The upward centroid movement was also augmented by the formation 
of a broad secondary layer above the main peak, as can be seen in the 
September and November observations (Figures 1 and 2). The occurrence 
of a minimum in mixing ratio between the main and secondary layers sug- 
gests that upward diffusion is not the sole source of the secondary 
layer, and that a particle-formation source may be present above the 
main layer. As noted in Section 3, Crutzen (1976) has shown that photo- 
dissociation of volcanically injected carbonyl sulfide (CSO) could pro- 
vide such a source. (See also the discussion in Section 4.c on possible 
normalization errors.) 

Compared to some previous estimates of the residence time of strato- 
spheric aerosols, the post-Fuego decay was rather rapid. For example, 
an "average residence time" or "average lifetime" of 1.5 to 2 years for 
stratospheric particles is frequently mentioned (e.g., Junge, 1974; 
Castleman, 1973), whereas the vertically integrated particulate back- 
scattering in 1975 decayed exponentially with a mean life of only about 
eight months. One reason for this apparent discrepancy is the relatively 
low altitude at which the Fuego injection occurred. Calculations by 
Hunten (1975) have recently emphasized the strong dependence of particle 
residence times on altitude, especially between the tropopause and abou^ 
20 km. (For the case modeled by Hunten, a decrease in altitude below 


(20 km by 3 km approximately halved the residence time. See also 
Reiter, 1975.) Further, Moore et al. (1973) concluded from measurements 
that stratospheric particles below 15 km have a residence time of only 
one to two months. Since the average 1975 centroid of particulate back- 
scattering was near 16 km, well below preinjection values (see Figure 4), 
it is not surprising that the 1975 residence time was small compared to 
values for a nominal "20-km aerosol." The short 1975 residence time was 
also partially caused by the size of the 1975 particles, which was, on 
average, larger than preinjection values (Ferry, 1975; Hofmann and Rosen, 
1976). 

For these reasons the present results are not necessarily at variance 
with previous estimates of the residence time of typical background 
stratospheric particles at altitudes near 20 km. They do, however, in- 
dicate that previous estimates of the duration of Influence of single 
volcanic injections may have been too large, at least for injections 
(e.g. , by Agung, see Figure 7) at altitudes appreciably below the unper- 
turbed aerosol peak. We note from Figure 5 that the decline in lidar 
backscattering measured in Australia after the post-Fernandina maximum 
had a time constant similar to that observed by both balloon and lidar 
after the post-Fuego maximum (also shown in Figure 5). Several other 
authors have already pointed out that the apparently much slower decline 
of the post-Agung maximum probably was the result of several subsequent 
smaller volcanic injections. 
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In the context of generalisations about the effects of volcanic in- 
jections, it should be noted that a change in dominant particle composi- 
tion, from sulfate to silicate, did occur after the Fuego eruption, 
although this has evidently been typical of previous eruptions (see, 
e.g. , Cadle and Grams, 1975; Pollack et al«, 1976a, b). (This fact has 
a bearing on our computations to follow.) Measurements made in the 
Fuego stratospheric plume very soon after the eruption revealed a particle 
composition that was predominantly sulfate, with no detectable silicate 
(Ferry and Lem, 1975). Four months after the eruption, a balloon measure- 
ment showed that essentially all of the stratospheric particles were 
volatile, again ruling out the possibility of a significant population 
of silicates (Hofmann and Rosen, 1976). Since the radiative properties 
of sulfates and silicates are significantly different, this lack of 
change in composition has a bearing on the radiative and climatic effects 
of the post-Fuego aerosol layer, as discussed below. 

b. Radiative and thermal effects 

The post-Fuego stratospheric aerosol layer scatters and absorbs both 
solar radiation and thermal infrared radiation; it therefore can modify 
radiation fields and temperatures in the atmosphere and at the surface. 
Measurements of the layer have determined particle concentration, com- 
position, and size distribution rather well, and several radiative, or 
radiative /convective, models are available to use these experimental 
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inputs in computing the radiative and climatic effects of the layer. Al- 
though such a computation is beyond the scope of the present study, some 
useful estimates of the radiative and climatic effects of the post-Fuego 
aerosol can be derived on the basis of several modeling results already 
published. While the particle characteristics assumed in the published 
models are not exactly those of the post-Fuego layer, the model results 
are nevertheless useful for gauging the possible size of effects and for 
guiding more detailed computations. 

Table 2 presents radiative and thermal properties and consequences 
of stratospheric aerosol layers, as computed by four recently published 
models. The detailed assumptions and methods used in each of the models 
are described in the references given in Table 2; a brief synopsis of 
these assumptions and methods is given in Table 3. To explain the symbols 
used in Table 2 we recount briefly the radiative and thermal impacts of 
a stratospheric aerosol layer, which is characterised in Table 2 by its 
mid-visible (X » 550 urn) optical thickness, t. 

We first consider the effects of the layer on solar radiation. At 
any given time and location a fraction R of incoming solar radiation is 
scattered by the layer back to space* The average of this fraction over 
all solar senith angles for all global locations is Rg, the global solar 
albedo of the stratspheric aerosol layer [row (1) of Table 2). The par- 
ticles and gases in the layer also absorb a fraction a of Incoming sun- 
light [row (2)]. The combined effect of reflection and absorption by 
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Tabla I 


MOUTHS MO TKMM CmCfl or tmtOiKIMIC AEROSOL Unit, 
AS MHBIC 1KB ST SEVERAL UCBCT MODELS 


Altec tad Par am tar 



Solar aftecta 


(0. )> t O.U)T 


(1) a (by particle.) , <0.02r 

a (by partialaa ' 
and (aaaa) 


: (0.1? t 0.07 ) At* 
| >(0.25 * 0.11) At* 


-0.2AAT* 


Tharaal infrared 
aftecta 


CoatbiMd aftecta 


M'l day )At (baaalt part tclea 
onlyl 


(-6.1'K to -10* K)i* 


Jt /(1 . D. 22 T)* IAS'S* 

( [<1.0.J?t)/(1 -0.11t)J. 1.9A « 10* U (t A ( K)]"*T ) 

“J l * I - 0 . 10 * j 

« IA2‘S* 


H» have auuMd an optically effective s»an particle radius of 0.3 t 0,2 ufl. 

♦ 

Sm also Toon and Pollack (1976). 

* 

Aaauaaa an unperturbed system albedo A * 0.3. and no effects of increased photon path laotth In absorbing gases. 
Computed Iron In, (3) and (A), 

* 

A* sines an unperturbed ayata* alba do A • 0.3. a no -part lei* T ■ 2U R« and R a, c as io rows (l), (5), and (6). 

A C 


ORIGINAL PAGE IS 
OP POOR QJAIJTY 









Table 3 


PARTKU CHARACTERISTICS AND MQHUK METHODS 
OBBS ffl THE HOBBLE CtlCD ff TKStS 2 



Hodal J 

Charactarlatlc 
or Matbod 

Cadla aad Oraaa 
(1979) 

Pollack at al. 
(1976)* 

Ptnnlok at al. 
(1976) 

Harehvardhao 
and Ceaa (1979) 

Coapoaldon 

Sulfuric acid* 

791 Hj 80 4 , 291 HjO 


791 Bj 80 4 , 291 HjO 

Rafracciva ladax 
(vlalbla) 

l.43-kl, k< 0.002* 

Palaar and Wllllana 
(1979) 

l.A).kl, k< 0.0027 

Palaar and Wllllana 
(1973) 

Rafracciva lndax 
(lnfrarad) 


PaUaar and Wllllana 
(1979) 


Palaar and Wllllana 
(1975) 

Sica dlatrlbutlon 

t 

told; r B • 0.0)5 ua. 
o s 2.0 

Lognormal; 
r. ■ 0.0725 ua. 
o > 1.36 

Baaabarg (1973) 

Vlalbla radlatlva 
cranafar 

Slngla acattarlng, 
partlclaa only 

Multiple acattarlng 
(Inhoaoganaoua dou- 
bling). partlclaa 

and gaaaa 

Slngla acattarlng. 
partlclaa only 

Multiple acattarlng 
(Coaklay and Cbylak 
(1973)1. partlclaa 
only 

lnfrarad radlatlva 
tranafar 

Pollack (1969). 
partlclaa and gaaaa 

Newtonian cooling coafficlanta, 

partlclaa and gaaaa 

Sorfaca taaparatura 
calculation 

Global anargy balanca, plua 
radlatlvfconvactlva equilibrium 
(few eaaaa) 

Clobal anargy balanca 


* 

Sat alto Toon and Pollack (1976). 

Salactad by ua from aavaral cholcaa given. 

Interpolated by ua from aavaral cholcaa glvtn. 

{ 

Raaulta praaentad for no nod la pa r parelclaa with a wide rang* of alaaa. 
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the layer is, in general, to change the albedo 4 of the aartfa-acroocp hag e 
aye tea by an amount AA Crow (3)3* Below die stratospheric layer, the 
effect of thla reflection and absorption la a loea of Incoming tolar 
radiation. In a eenee thla loaa may be conaldered aa equivalent to a 
reduction of the aolar conatant by a fraction 18/1 (row <4)1. 

Since the 8tratoapherlc aerosol layer la optically thin (t « l), 

the models yield that its albedo Bg and ita fractional absorption a are 

* 

directly proportional to its (mid-visible) optical thickness t [See 
rows (1) and (2) of Table 2]. The system albedo change AA is, in the 
more realistic models, computed from the equation of radiative transfer 
by comparing the ratio of incoming and outgoing solar lrradiance for 
perturbed and unperturbed systems. In the simplest models, which con- 
sider an isolated aerosol layer in an otherwise transparent atmosphere, 
this albedo change may be simply computed as (e.g. , Bussell and Grams, 
1975): 

AA - A' - A - a£(1 - Bg - a) 2 /(l - ARg) j + Rg - A . (3) 

In either case, the albedo change again turns out to be directly propor- 
tional to layer optical thickness t, or to the change At with respect 
to a background or unperturbed layer (See row (3) of Table 2). Similar 
considerations apply to the effective reduction in solar constant AS/S, 
which, for the simplest models, Is given by (e.g., Cadle and Grams, 1975): 
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A 8/8 « -Rgd - A) - a . (4) 

Again, for all models, the change is proper clonal to t or At (goo row (4) 
of Table 2). 

The aerosol layer also has affects in the the nasi infrared, where 
it both absorbs and emits energy (scattering is practically negligible)* 
Only In the model of Harshvardhaa and Cess are the absorption ft and emit* 
tance « explicitly computed; they also are directly proportional to 
t [rows (5) and (4) of Table 2]. In the model of Pollack et al. (1976a) 
infrared effects are determined by a solution of the radiative transfer 
equation (Pollack, 1969)* Infrared effects of sulfuric acid particles 
are not considered in the models of Cadle and Crams (1975) and of Pin- 
nick et al* (1976), 

The combined effects of the stratospheric aerosol layer on solar 
and infrared radiation in general alter temperatures both in the vicinity 
of the aerosol layer itself and at the earth's surface. Changes at the 
aerosol layer are most pronounced in the months near the postvolcanic 
peak of aerosol concentration, because of the relatively short time con- 
stant (— 100 days) of the free atmosphere for temperature changes 
(Pollack et al. , 1976a, b). Changes at the earth's surface take longer 
to develop because of the long time constant (— 4 years) of the earth- 
ocean system. Pollack et al. (1976a) compute both perturbed heating 
rates AH A at the aerosol layer and surface temperature changes AT g . For 


32 


optically thin layer* both changes torn out to bo proportional to tho 
aid-visible optical thickness change At (gee rows (?) sad (9) of Table 21* 
(Unfortunately, because many paet volcanic particulate injections have 
initially consisted predominantly of silicates, the AB A computation was 
done for basalt particles only.) Barshvardhan and Ceos (1975) provide 
explicit expressions for the change* AT A , AT g in both aerosol and surface 
temperatures as a function of t [ See row* (8) and (9) of Table 2]. Al- 

i 

though formally nonlinear, these expressions also turn out to be practi- 
cally linear for the optical thickness changes At with which we will deal. 

To explore the radiative and thermal consequences of the post-Fuego 
stratospheric aerosol layer predicted by the four models of Table 3 we 
adopt a odd-visible optical thickness t of 0.030, tdiich is the nine- 
month mean of our 1975 observations (Table 1). We also adopt a non vol- 
canic, or unperturbed, optical thickness of 0.005, which Is the unper- 
turbed value assumed by Pollack ct al. (1976a) and which falls within 
the range of both our pre-Fuego determinations and those of Pinnlck et 
al. (1976) (Table l). The results of substituting these values into the 
expressions of Table 2 are shown In Table 4. 

The three analytical models presented in the table (Cadle and Crams 
(1975); Pinnlck et al. U976); Harshvardhan and Cess (1975)1 agree in 
predicting a layer albedo Rg of about 1%; all four models show that 
solar absorption by the layer is negligible. The three analytical models 
predict an Increase AA in system slbedo that is about 1.5% of the 
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TabU A 


RADIATIVE AND THERMAL EFFECTS OF A STRATOSPHERIC 
AEROSOL LAYER WITH MtD-VXSIBLB OPTICAL THICKNESS 
T - 0.030 AS PREDICTED BY THE MODELS IN TABLES 2 AND 3* 



Aaautaea an unperturbed mld-vlalble atratoapherlc aeroaol optical thlckneaa * » 0.005; 
hence Lr « 0.025. Changea in rowa (2), (4), (7), (8). (9) an with ratpect \z the unper- 
turbed ayaten. 

^We have aaauned an optica I '.v effective nv’n particle radlua of 0.3 : 0.2 „». 

% 

See also Toon and Pollack (1967). 

r Aeotmea an unperturbed avatem albedo A « 0.3, and no effecta of increaaed photon path 
length in absorbing gaaea. 'Jomputed from Eqa. (3) and (4). 

** 

Aaauraea an unperturbed ay a tern albedo A * 0,3, a no-particle T * 214 K, and R., a, and < 
aa in rowa (1). (5), and (6). 0 
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unperturbed albedo (A ■ 0.3). The prediction of Pollack at al. (1976a) 
is somewhat less than half of this value, possibly because their model 
accounts for increased solar absorption caused by longer photon path 
lengths in the aosorbing gases of the aerosol layer. There is a similar 
difference between the numerical and analytical models in the predicted 
change of sublayer insolation. The analytical models predict a decrease 
of about 0.7%, the- americal model only about 0.35%. 

Perhaps of more interest are the predicted changes of atmospheric 
and surface temperatures caused by the combination of solar and infrared 
effects. For an optical thickness increase of 0.025, the model of Pollack 
et al. (1976a) predicts an increased heating rate AH^ at the altitude 
of the layer, of about 0.05? K per day. Such an effect, if persistent 
for several months, would produce observable temperature increases. Un- 
fortunately, the Pollack et al. (1976a) calculation was made only for 
basalt particles, and so does not apply specifically to the post-Fuego 
particles, which were sulfuric acid. On the other hand, the calculations 
of Harshvardhan and Cess (1975) are for sulfuric acid particles. Their 
model predicts an equilibrium aerosol temperature increase AT^ of about 
2.4°K, which again might be observable by comparing average temperatures 
for the year before and after the Fuego eruption, or other suitable 
periods. The fact that the first three or four months of 1975 had op- 
tical thickness values significantly in excess of our 1975 mean of 0.03 
makes such a search for a heating effect appear even more worthwhile. 
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(Note, however, the reservation on particle size distribution mentioned 1 
below.) Preliminary observations of such a heating effect have in fact 
been reported by McCormick etal. (1973) and Hirono et al. (1973). While 
these preliminary indications appear to be in qualitative accord with 
the model predictions, we do caution that stratospheric temperature changes 
caused by dynamic and other effects (especially breakdowns in the winter 
polar vortex) could overwhelm or appear to augment the purely radiative 
changes predicted here. See, for example, Belmont (197S). 

Both the models of Pollack et al. and of Harshvardhan and Cess 
predict a negative change ATg in equilibrium surface temperature for in- 
creased stratospheric aerosol optical thickness. (That is, the models 
predict that the increased solar albedo effect dominates over the in- 
creased infrared trapping or "greenhouse" effect, although the infrared 
effect is significant.) For a persistent layer of optical thickness 
0.030 the Harshvardhan and Cess model predicts a surface temperature de- 
crease of 0.8°K, whereas the Pollack et al. model prediction is about a 
factor of four less. Although a decrease of 0.8°K in surface temperature 
is on the threshold of detectability, it must be recalled that the long 
response time < •— 4 years) of the ocean-earth system will substantially 
diminish this decrease in the first year, and that the rapid removal rate 
of the post-Fuego aerosol burden (t ** 8 months) will rapidly eliminate 
the decrease thereafter. Feedback mechanisms (principally extension of 
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snow and Ice cover) chat are not considered in either the Harshvardhan and 
Cess or Pollack et a*« models, are also relatively slow, so that their am- 
plification o£ radiative effects, as predicted hy some other models 
(e.g., Budyko, 1969; Sellers. 1973) would evidently not be operative in 
the Fuego case. Thus it seems unlikely that a demonstrable connection 
between any recent surface temperature changes and the post-Fuego strato- 
spheric aerosol could plausibly be made. This result, of course, in no 
way rules out the general possibility of volcanic influence on surface 
climate; Pollack et al. (1976b) have shown that .he results of larger 
eruptions (e.g.. Krakatoa, Katmai, Agung), or of more frequent ones, 
could indeed be climatically significant. 

While we feel that the results presented above are instructive, we 
also point out that the calculations can and should be improved. For 
example, the measured post-Fuego size distribution (Ferry, 1975; Hofmann 
and Rosen, 1976) should be used in place of the distributions listed in 
Table 3. It would also be instructive to repeat the Pollack et al. com- 
putation of AH^ using a sulfuric acid (rather than basalt) composition 
model. And, finally, since volcanic injection episodes are decidedly 
time -dependent, time -dependent computations for AT^ and ATg, using the 
appropriate system response times, would be useful. Of course, the 
reasons for differing model predictions should also be explored and re- 
solved. 
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c. Accuracy of the Ildar results 


As described in Section 3, all of our particulate backseat taring 
values are derived by normalizing each scattering ratio profile to force 
its minimum value to unity. This procedure yields valid profiles pro- 
vided an altitude with negligible ratio of particulate to gaseous back- 
scattering exists within the altitude range over which signals are ana- 
lyzed. If such a level does not occur, then particulate backscattering 
is underestimated at all altitudes (and this error is not included in the 
error bars of Figures 1 and 2). The magnitude underestimation 

cannot be determined from the lidar data alone, but it can be estimated 
by use of independent measurements of the vertical profile of particle 
number mixing ratio, as derived, for example, from dustsonde measurements 
(Hofmann et al., 1975, Hofmann and Rosen, 1976). We have previously 
used pre-Fuego dustsonde measurements to demonstrate that the average 
error introduced into our pre-Fuego backscattering values by the normaliza- 
tion procedure was less than 10% of the peak backscattering measured in 
a typical profile (Russell et al, 1976b). A similar inspection of post- 
Fuego dustsonde mixing ratios (Hofmann and Rosen, 1976) indicates that 
the average normalization error for our post-Fuego backscattering mea- 
surements was also less than 10% of profile peak backscattering values. 

It is true that the tropopause region in 1975 was frequently quite dirty, 
forcing the normalization to altitudes above the main peak, where optical 
mixing ratio if frequently nonnegligible; nevertheless, the very large 
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Increase in peak backseat taring during the post-Fuego period prevented 
the ratio of error to peak backseat taring from increasing in the post* 
Fuego period. Thus we feel that to compensate for normalisation errors 
any increase in the post-Fuego optical thickness values of Table 1 must 
be less than 10% of the listed values. 

In connection with the discussion of a possible secondary layer 
above the main peak, it should be pointed out that the normalisation 
procedure could also lead to an apparent, but artificial, increase of 
scattering above the main peak when scattering below the main peak de- 
creases, as it did during 1975. This follows because, when scattering 
ratios below the main peak are larger than anywhere above the main peak-- 
as was typical in early 1975 — the normalization procedure forces the 
scattering ratio to unity (implying no particulate scattering) at some 
point above the main peak, even though actual particulate scattering may 
be nonzero there. Then later, as low-altitude scattering ratios decrease, 
the normalization altitude would move below the main peak, and scattering 
above the main peak could appear to increase while actually remaining at 
a constant, nonzero value. While we feel that this possibility is a real 
one and must always be kept in mind, we point out that the shapes of 
scattering ratio profiles are completely independent of the normaliza- 
tion procedure. Our statements regarding the development of a secondary 
layer are based primarily on the shapes of scattering ratio profiles-- 
speciflcally the minimum just above the main peak and the broad, but 
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well-defined, hump above it, as shown in Figure 2. We therefore feel 
that these statements are valid regardless of any normalisation ambigu- 
ities. 

We have also noted in Section 2 that the Elterman (1968) extinction 
model was assumed in computing 1975 scattering ratio profiles. This 
model is in fact listed in the first column of Table 1. Comparison with 
the last column of Table 1 shows that, at any altitude above 10 km, the 
difference in two-way, 700-nm optical thickness between the model and 
the 1975 mean measurements was less than 0.01 (again assuming extinction 
a \ *). This difference is the maximum relative change in scattering 
ratio that would result from using the last column of Table 1 in place 
of the first and recomputing. Even for the very last measurements made 
in 1975, when the difference between model and measurements was about 
twice this large, the change would have been less than the error bars 
shown in Figures 1 and 2, which are about ± 3%. The resulting change in 
vertically integrated backscattering (hence turbidity and optical thick- 
ness) would have been negligible. 
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5* Summary 

Our 1975 stratospheric Ildar observations began at about the peak 
of northern midlatitude influence of the Fuego volcanic injection. 

Multiple- layered structure, observed by others since the time of the 
eruption, was still present in early 1975, but by May the individual 
layers had merged to form a single, broader layer. The increase in par- 
ticulate backscattering resulting from the volcanic injection was confined 
almost entirely to altitudes below 20 km. As a result, layer centroids 
in 1975 were typically several km below their altitudes prior to the 
eruption. After April 1975 the centroids gradually rose toward their 
preinjection value, as particles at the lowest stratospheric altitudes 
were removed faster than those at higher altitudes. Toward the end of 
our 1975 observations, profile shapes suggest the formation of a secondary, 
broad layer above the main peak. A search for dawn-to-twilight differ- 
ences in layer structure, and for differences caused by the rising sun, 
revealed no changes that were systematically different from those commonly 
observed during other time periods of similar length. 

From late February on, both vertically integrated particulate back- 
scattering and the peak ratio of particulate to gaseous backscattering 
displayed approximately exponential declines, with mean lifetimes of 
8 and 11 months, respectively. Very similar behavior was observed by 
balloon-borne particle counters in Wyoming during the same period 
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(Hofmann and Rosen, 1976). The relatively abort residence tines are a 
combined result of the low altitude of the voleanlc injection and of the 
larger mean particle siae of the pose volcanic aerosol with respect to 
the prevolcanic or "background" aerosol. 

The 1975 mean profile of scattering ratio had a maximum value of 

1.7, and vertically integrated particulate backseat ter lng (10*30 km) of 
*4 *1 

3.6 x 10 sr (both for X « 694 nm). When combined with several real* 
Istlc optical models of the sulphuric acid particles, this backscatterlng 
yields a mean mid-visible (X « 550 tun) optical thickness of 0.03 ± 0.005. 
This value is about a factor of six larger than the preinjection optical 
thickness, but smaller than values observed in northern midlatitudes for 
several years following the Agung eruption of 1963. 

Radiative and thermal consequences of the measured post-Fuego layer 
were computed using several recently published models. The models yield 
a global albedo for the 1975 mean layer of about 1%, or possibly smaller 
because of gaseous absorption of photons scattered in the layer. The 
derived increase in earth -atmosphere system albedo was between 0.002 
and 0.005. The models predict that infrared absorption by the sulfuric 
acid particles could Increase temperatures at the aerosol layer altitudes 
by several degrees K. They predict a decrease in surface temperature 
(that is, increased solar albedo dominates over increased infrared trap- 
ping), but the decrease is considerably smaller than 1°K, both because of 
the small optical thickness of the volcanic layer and because of its short 
residence time relative to the earth-ocean thermal response time. 
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